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In recent experiments with ion traps, long-range interactions were associated with the exceptionally fast 
propagation of perturbation, while in some theoretical works they have also been related with the suppression of 
propagation. Here, we show that such apparently contradictory behavior is caused by a general property of long- 
range interacting systems, which we name Cooperative Shielding. It refers to shielded subspaces that emerge as 
the system size increases and inside of which the evolution is unaffected by long-range interactions for a long 
time. As a result, the dynamics strongly depends on the initial state: if it belongs to a shielded subspace, the 
spreading of perturbation satisfies the Lieb-Robinson bound and may even be suppressed, while for initial states 
with components in various subspaces, the propagation may be quasi-instantaneous. We establish an analogy 
between the shielding effect and the onset of quantum Zeno subspaces. The derived effective Zeno Hamiltonian 
successfully describes the short-ranged dynamics inside the subspaces up to a time scale that increases with 
system size. Cooperative Shielding can be tested in current experiments with trapped ions. 

PACS numbers: 03.65.Xp; 75.10.Pq; 37.10.Ty; 67.85.-d 


Introduction.- A better understanding of the nonequilib¬ 
rium dynamics of many-body quantum systems is central to a 
wide range of fields, from atomic, molecular, and condensed 
matter physics to quantum information and cosmology. New 
insights into the subject have been obtained thanks to the re¬ 
markable level of controllability and isolation of experiments 
with optical lattices M and trapped ions | f| [§j]. Recently 
there has been a surge of interest in the dynamics of sys¬ 
tems with long-range interactions, triggered by experiments 
with ion traps |@, 9], where the range of interactions in one¬ 
dimensional (ID) spin models can be tuned with great accu¬ 
racy. Other realistic systems that contain long-range interac¬ 
tion include cold atomic clouds 0, natural light-harvesting 
complexes ITTI- L3ll . helium Rydberg atoms 0, and cold Ry¬ 
dberg gases llal . Long-range interacting systems display fea¬ 
tures that are not often observed in other systems, such as 
broken ergodicity EMI and long-lasting out-of-equilibrium 
regimes [20]. 

According to the usual definition f2lll , in d dimension, an 
interaction decaying as l/r a (where r is the distance be¬ 
tween two bodies), is short range when a > d and it is long- 
range when a < d. A major topic of investigation has been 
whether the propagation of excitations in systems with long- 
range interaction remains or not confined to an effective light 
cone ll22l - [30ll . as defined by the Lieb-Robinson bound IpT I] 
and its generalizations (|3().| and references therein). In the 
aforementioned experiments with trapped ions, it was ob¬ 
served that for short-range interaction, the propagation of per¬ 
turbation is characterized by a constant maximal velocity, 
being bounded to an effective light cone. As a decreases, 
the propagation velocity increases and eventually diverges. 
For long-range interaction, a < 1, the light-cone picture is 
no longer valid and the dynamics becomes nonlocal. How¬ 
ever, examples of constraint dynamics in long-range interact¬ 


ing systems have also been reported, including logarithmic 
growth of entanglement |23], light-cone feature s_|130ll . self¬ 
trapping i32fl , and slow decays at critical points 13311 . 

Here, we show that these contradictory results are due to 
a general effect present in long-range interacting systems, 
which we name Cooperative Shielding. It corresponds to the 
onset of approximate superselection rules that cause a strong 
dependence of the dynamics on the initial state. Inside a su¬ 
perselection subspace, long-range interactions do not affect 
the system evolution ( shielding ) up to a time scale that grows 
with system size ( cooperativity ). The dynamics can then be 
described by an effective short-ranged Hamiltonian that either 
leads to a propagation within the Lieb-Robinson light cone or 
to localization. In contrast, for an initial state with compo¬ 
nents over several subspaces, the propagation of excitations is 
affected by long-range interactions and can be unbounded. 

To explain how shielding can arise in a very trivial case, 
let us consider the total Hamiltonian H = Hq + V, de¬ 
scribing a many-body quantum system, where Hq has one- 
body terms and possible short-range interactions, and V cor¬ 
responds to some additional interactions. If [Hq, V] = 0 and 
V is highly degenerate in one of its eigensubspaces V, so 
that V\V k ) = v\Vk) V|ufc) 6 V, the evolution of any initial 
state |^o) belonging to such eigensubspace is simply given 
by: |^(f)) = e~’ lvt l h e~ lHat l h \tl}Q). Since the only effect 
of V is to induce a global phase, the dynamics is shielded 
from V and determined only by Hq. In contrast, if the initial 
state has large components in more than one eigensubspace of 
V, the dynamics will not be shielded from V. The question 
that we now pose is whether shielding is still possible when 
[Hq, V] f 0 and V is no longer degenerate. We show that the 
answer is positive when V involves only long-range interac¬ 
tions. The dynamics can remain shielded, but now for a finite 
time that increases with system size. 
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One can also draw a parallel between the picture above and 
the quantum Zeno effect (QZE). In the QZE, the dynamics 
of the system remains confined to subspaces tailored by the 
interaction with a measuring apparatus 1341 - 13811 . The stronger 
the interaction is, the better defined the subspaces become. 
Here, instead, the interaction strength is kept fixed, but due to 
its long-range-nature, invariant subspaces are generated. The 
dynamics, restricted to the invariant subspaces, is described 
by a short-ranged Zeno Hamiltonian up to a time scale that 
diverges with system size. 

The Model.-We consider a ID spin-1/2 model with L sites 
and open boundary conditions described by the Hamiltonian, 


H = H 0 + V, 

L L — l 

H 0 = + h n )a ~ ^ J z a~ n (Jn+ii 


( 1 ) 




J 
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spin polarization, (cr“(t)), for an initial state where all spins 
point up in x, except for the spin in the middle of the chain, 
which points down, so M x = L/2 — 1 and 6=1. 
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Above, h = 1 and a^ ,y,z are the Pauli matrices on site n. The 
transverse field has a constant component B and a random part 
given by h n , where h n £ \—W/2,W/2] are random numbers 
from a uniform distribution. The nearest-neighbor (NN) in¬ 
teraction in the z-direction, of strength J z > 0, may or not be 
present. J is the strength of the interaction in the x-direction 
with a determining the range of the coupling. Unless spec¬ 
ified otherwise, J = 1. The Hamiltonian with W = 0 and 
J z = 0 describes the systems studied with ion traps |8. [f|. In 
agreement with those experiments, where a limited range of 
system sizes is explored, V is not rescaled by L. 

When a = 0, H can be written in terms of the total x- 
magnetization, M x = J2n-i °n/2, as 

L L — l jjr 

H = ]T(B + + £ J*«+1 + 2 JMl - —. (2) 

n=1 n=1 

The spectrum of V is divided into energy bands, each one 
associated with a value of the collective quantity M x . Each 
band, with energy Ef, = 2J(L/2 — 6) 2 — JX/2,has states with 
b and L — b excitations, where 6 = 0,1,... L/2. For instance, 
6=1 corresponds to states with one spin pointing up in the 
x-direction in a background of down-spins or vice-versa. An 
energy band contains 2 ([') degenerate states if 6 < L/2 and 
(k) states when 6 = L/2. In contrast, for 0 < a < 1, the 
states in each band V are not all degenerate anymore. 

Light Cones.- In Refs. Hi , the acceleration of the spread¬ 
ing of excitations and eventual surpassing of the Lieb-Robison 
bound achieved by decreasing a was verified for initial states 
corresponding to eigenstates of // (J , where each site had a spin 
either pointing up or down in the z-direction. These initial 
states have components in all subspaces of V. 

Motivated by the special role of the x-direction in Eq. Q 
and to show the main features of Cooperative Shielding, here 
we change the focus of attention to initial states with spins 
aligned along the x-axis. They are the eigenstates of V and 
are denoted by |I4). In Fig.Q] we show the evolution of the 


FIG. 1: (Color online) Density plots for the evolution of (o 
L = 13; B = 1/2; W = 0. Initial state: ( 0 - 7 ( 0 )) = —1 and 
(o-^j 7 (0)) = +1. A light cone typical of short-range interaction is 
seen in (a), as expected, but also in (d), (e), and (f) where the evo¬ 
lution is shielded from the present long-range interaction. Freezing 
occurs for very long times in (b); it also happens in (c) where the 
bands of V are not degenerate. 

In Fig. [I] (a), where the interaction is short range (a = 3), 
Ho effectively couples states belonging to different subspaces 
of V. The effects of both H 0 and V lead to the evident light 
cone. This is no longer the case for long-range interaction 
(i a < 1), as exemplified in Figs.[T](b) and (c) for a = 0 and 
0.5. Their dynamics is frozen for a long time, which increases 
with the range of the interaction [compare the time scales in 
(b) and (c)] and with the system size (see discussion below). 
The long-time localization of spin excitations in Figs. □ (b) 
and (c) is caused by both combined factors; the separated en¬ 
ergy bands of V and the absence of direct coupling within the 
band (Hq is not effective and J z = 0). Notice that the energy 
bands for case (c) are no longer degenerate, yet localization 
persists for a long time. 

Since the initial state is not an eigenstate of the total Hamil¬ 
tonian, the spin excitation does eventually spread and the spins 
reverse their signs (see Figs.^(b,c) and discussion in 14111). 
This magnetic reversal can be explained in terms of macro¬ 
scopic quantum tunneling lfl7h . 

While for a < 1 in the presence of an external field the 
dynamics is frozen, the addition of NN interaction ( J z ^ 0) 
restores the propagation of perturbations [Figs. □(d), (e), (f)]. 
Despite the existence of long-range interactions, the evolu¬ 
tion can be described by an effective short-ranged Hamilto¬ 
nian, as we show below. This is the hallmark of the Cooper¬ 
ative Shielding effect discussed in this work, the suppression 
of propagation [Figs, □(b), (c)] being only a special case of it. 

In Figs. □ (d), (e), (f), a light cone typical of short-range 
interactions emerges: the dynamics is independent of system 
size and of the long-range coupling J. In Fig. □(!'), J is twice 

















3 


as large as in Figs. |T|(d,e), but the results in the three panels 
are very similar, apart from border effects. The propagation of 
excitations depends only on J z up to long times. This shielded 
evolution occurs for any a < 1 (see more figures in I4dl ). In 
the case of 0 < a < 1, as in Figs.[T|(e) and (f), the bands of 
V are no longer degenerate, so the various eigenstates of V 
that are excited within the band have different eigenenergies. 
One could then expect V to affect the evolution, yet the veloc¬ 
ity of propagation remains independent of V for long times. 
This shows that the cause for shielding is not only the suppres¬ 
sion of the transitions between different bands of V, but also 
the narrow distribution of the energies of V inside the band. 
The motion remains constrained to subspaces that are quasi¬ 
degenerate w.r.t. to V. The emergence of quasi-constants of 
motion is recurrent in long-range interacting systems 0. 

Invariant Subspaces and Zeno effect.- Stimulated by the 
results of Fig. [7] we now analyze in more details the effects 
of infinite-range interaction (ct = 0) and their dependence on 
system size. For a general treatment, we assume a random 
transverse field, so B = 0 and h n ^ 0. We take as initial state 
|\I>(0)) a random superposition of all states \V^) that belong 
to the same fixed band b chosen for the analysis. We verified 
that the results for single states \V£) picked at random from 
the same energy band are equivalent. 

In Figs. [2] (a) and (b), we compute the probability, 
for the initial state to remain in its original energy band b, 

n(i) = EK^I e “ iif ^(0))l 2 ’ (3) 

k 

where the sum includes all the states of the selected energy 
band. The results are shown for ( Pt(t )}, where (.} indicates 
average over random realizations and initial states. We show 
the case of b = 1, but similar results hold for other bands. 
It is evident that the probability to remain in the initial band 
increases with system size. This happens in the presence of a 
random transverse field [Fig. [2] (a)] and also when NN inter¬ 
actions are added [Fig. 13(b)]. 

In Figs. 0(c) and (d), we plot the asymptotic values of the 
leakage probability, Pjeafc = l~lim t _ s . 00 (Pb(f)), as a function 
of the random field strength for J z = 0 [Fig. [2](c)] and v,v the 
NN coupling strength for W = 0 [Fig.[2](d)]. Pi ea k represents 
the probability for |\I/(0)) to leak outside its original band. It 
decreases with L, showing that as the system size increases, 
the evolution of | (0)) remains more and more confined to a 

subspace of V for a longer time. Note that the distance be¬ 
tween the bands nearby the initial one increases with L, but 
so does the number of states which are connected by Hq. The 
suppression of leakage takes into account this non-trivial in¬ 
terplay. A perturbative argument leads to P[ ea k ex (W/ J) 2 /L 
for W ^ 0 and J z = 0, while Pi ea k ex (J~/J) 2 /L forNN in¬ 
teraction only 61- Such scaling relations are consistent with 
our numerical data in Figs.[2](c) and (d). 

The invariant subspaces generated by long-range interac¬ 
tion can be related to the QZE ll34l - l38ll . This term refers to the 
familiar freezing of the dynamics due to frequent measure¬ 
ments, but also to the onset of invariant Zeno subspaces that 
occurs in unitary dynamics due to strong interactions [i36ll38?l 
and which has been studied experimentally |42|]. The latter is 
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FIG. 2: (Color online) Probability for the initial state to remain in 
(a,b) or leave (c,d) its original energy band. In (a,b): (Pj,(t)} for 
the initial random superposition of states \V£) from band 6 = 1 for 
L = 10,12,14 from bottom to top; in (a): J z = 0, W = 2 and in 
(b): J z = 1, W = 0. In (c,d): Pi ea k vs W for J z = 0 (c), and v,v J z 
for W = 0 (d). Symbols represent numerical results and full lines, 
analytical estimates 62] with an overall fitting multiplicative factor. 
In all panels: averages over 50 realizations, B = 0, a = 0. 


closer to our case and can be explained as follows. Consider 
the total Hamiltonian H = H s + gH meas , which one may 
interpret as a quantum system described by H s that is contin¬ 
uously observed by an “apparatus” characterized by gH meas . 
In the limit of strong coupling, g —> oo, a superselection rule 
is induced that splits the Hilbert space into the eigensubspaces 
of J/meas- Each one of these invariant quantum Zeno sub¬ 
spaces is specified by an eigenvalue Vk and is formed by the 
corresponding set of degenerate eigenstates of H meas . The 
dynamics becomes confined to these subspaces and dictated 
by the Zeno Hamiltonian Hz = Yk HkH s Ilk + v k 11 a , where 
n fe are the projectors onto the eigensubspaces of iT meas cor¬ 
responding to the eigenvalues Vk- 

For the system investigated here, we associate H s with // () 
and .giJmeas with V. The subspaces of V, with fixed num¬ 
bers b of excitations, become invariant subspaces of the to¬ 
tal Hamiltonian not only when J —> oo with B. IT, J z fixed, 
which is the scenario of the QZE described above, but also in 
the large system size limit, L —> oo, which is the main focus 
of this work. 

When J z = 0, the Zeno Hamiltonian coincides with 
V, because the transverse field does not couple directly 
states \V k ) that belong to the same eigensubspaces of V, so 
Y k n fe 7J 0 n fc = 0. This explains why the dynamics in Fig. |T] 
(b) is frozen for very long times. On the other hand, in the case 
where B, W = 0 and J z ^ 0, we can rewrite Hq in terms of 
the cr „ x operators that flip the spins in the ^-direction. The 
projection of the NN part of the Hamiltonian on the eigen¬ 
subspaces of V leaves only the term Un x<J n+i + a n X( 7 n+ 1 > 
which leads to a Zeno Hamiltonian with an effective NN in¬ 
teraction that conserves the number of excitations inside each 
band b. This explains why in Fig. |T] (d) a light cone typical of 
short-range interactions appears. 

Fidelity Decay- To substantiate that the dynamics in the 
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subspaces with fixed b becomes indeed controlled by the Zeno 
Hamiltonian as L increases, we analyze the fidelity between 
an initial state evolved under the total Hamiltonian H and the 
same state evolved under Hz, 

F(t) = \{^{Q)\e iHzt e- iHt \^{Q))\ 1 2 . (4) 

It is clear that if H —> Hz then F(t) —> 1. The results are 
shown in Fig. 0 Equivalently to Fig. |2] we fix B = 0 and 
deal with averages over disorder and initial states, which gives 
( F(t )). |(0)) is again a random superposition of all states 
|I4 b ) belonging to the same band b. 



FIG. 3: (Color online) Fidelity decay and time for it to reach the 
value 1/2; initial states are random superpositions of |Vj?). Upper 
panels: F(t) for b = 3 for J z = 0, W = 2 (a) and for W = 0, J z = 
1 (d). From bottom to top: L = 10,12,14. Numerical results: full 
lines. Gaussian decay: dashed lines. Lower panels have J z = 0 and 
give T\j 2 Vi W for L = 12 (b), and vs L for W = 2 (c), for |\P(0)) 
from different bands. Numerical data: symbols. Analytical estimate 
Ti /2 = Ci/SE with ci a fitting parameter: dashed lines. All panels: 
averages over 50 realizations, a = 0, B = 0. 

In Figs. [3] (a) and (d) the fidelity is plotted vs time for dif¬ 
ferent system sizes for the band with b = 3. In panel (a), Hq 
contains only the random fields, while in (d), Hq contains only 
NN interaction. In both cases the fidelity decay slows down 
as the system size increases, confirming that Hz determines 
the dynamics for large L. 

For the J z = 0 case of Fig. [3] (a), since the projection of 
//(i on the b subspace is zero, the fidelity coincides with the 
survival probability, F(t) = |('I , (0)|'T(i))| 2 , which, counter¬ 
intuitively, decays slower as the system size increases. This 
shows that the dynamics localizes as L —> oc. F(t) decays as 
a Gaussian 143144711 - see dashed lines in Fig. [3](a). 

In Figs. [3](b) and (c) we study how the time Tj/ 2 that it 


takes for the survival probability to reach the value 1/2 de¬ 
pends on the disorder strength (b) and on system size (c). Fig- 
ure[3](b) provides information associated with the usual QZE, 
where the quantum Zeno subspaces are induced by decreasing 
the strength of Hq- One sees that the dynamics slows down 
with the reduction of disorder as (Tj/ 2 ) oc W~ 2 . In Fig [3] 
(c), (Tj/ 2 ) grows with L, corroborating our claims that the fi¬ 
delity increases and the excitations become more localized as 
the system size increases. 

The estimation of the dependence of Tj/ 2 on the param¬ 
eters of H goes as follows. Since the eigenstates of V in 
each invariant subspace are degenerate, the perturbation Hq 
mixes them all. In this case, the energy uncertainty u> of the 
initial state can be approximated by the energy spread 8E of 
each band induced by the perturbation. The fidelity decay can 
then be estimated as Tj/ 2 — 1/SE, where SE is computed 
from perturbation theory ||4T)| . For large system sizes one 
has Tj/2 oc Jy^L/W 2 . The analytical estimates for Tj / 2 are 
shown with dashed curves in Figs. 0(b) and (c). The agree¬ 
ment is excellent. 

We note that Tj / 2 gives the time scale over which the 
shielding effect persists. In finite systems, shielding is effec¬ 
tive for a finite time that can, however, be exceedingly long, 
as shown in Fig. [3] 

Conclusions.- We revealed a generic effect of long-range 
interacting systems: Cooperative Shielding. It refers to invari¬ 
ant subspaces that emerge as the system size increases. Inside 
these subspaces, the dynamics occurs as if long-range inter¬ 
action was absent, being dictated by effective short-ranged 
Hamiltonians. A parallel was established between these 
Hamiltonians and Zeno Hamiltonians. 

The analysis and control of nonequilibrium dynamics can 
never be detached from the initial state considered. For ex¬ 
actly the same Hamiltonian with long-range interaction, an 
initial state with components in the various subspaces induced 
by that interaction leads to a nonlocal propagation of pertur¬ 
bation, as demonstrated experimentally with ion traps SI, 
while an initial state belonging to a single subspace is unaf¬ 
fected by the long-range interaction, as verified here. Cooper¬ 
ative Shielding could also be tested by those experiments. 
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I. INTRODUCTION TO THE SUPPLEMENTARY 
MATERIAL 

Here, we provide further illustrations, reinforcing that 
shielding is a generic property of long-range interacting sys¬ 


tems. We also analyze the magnetic reversal of the spins in 
the presence of long-range interactions and show how we ob¬ 
tained our estimates for the leakage probability and the en¬ 
ergy spread SE of each band. The energy spread is used for 
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approximating T\/ 2 — 1/ SE, which refers to time for the sur¬ 
vival probability to reach the value 1/2. 


II. SHIELDING EFFECT 


To further support that in the presence of long-range inter¬ 
action (a < 1), the dynamics is shielded from V and therefore 
does not depend on the long-range coupling strength J for a 
long time, we show in Fig. [4] the evolution of the spin polar¬ 
ization in the ^-direction for different values of J. The initial 
state has the spins aligned in the x-direction as in the Fig.l 
of the main text. We fix the strength of the NN interaction to 
J z — 1. The figure shows that the speed of the propagation 
remains unchanged, as J increases from top to bottom. It de¬ 
pends only on the strength J z of the NN coupling. The same 
behavior is observed also as a increases from left to right. 


(a) J=OJ, a=0 jO (b) J=OJ, a=0J (c)J=OJ,a=OJ 



(d) J=UQ, a=00 it) J=t O. a =03 (f) J=I O. a=03 



(g) J=2.0, a=0.0 (h) J=2.0, a =03 (i) J=2.0. a =05 



1 3 5 7 9 11 13 1 3 5 7 9 11 13 1 3 5 7 9 11 13 

Site Site Site 


FIG. 4: (Color online) Evolution of the polarization (a%(t)) for all 
sites n; L — 13, B = 1/2, W = 0, and J z = 1. J increases from 
top to bottom and a from left to right. Initial state: (crf i+1 w 2 (0)) = 
— 1 and {&n^( L +i)/ 2 (0)) = Tl- 


HE MAGNETIC REVERSAL 


The frozen dynamics seen in Figs.l (b) and (c ) of the main 
text, for S / 0 and J z = 1, holds for a finite time. Eventually 
the spins reverse their signs. The polarized-sign reversal for 
a < 1 is a collective effect: the external field rotates all spins 
synchronously in the xy-plane, as seen in Figs. 0(a) and (b) 
below. 

In Fig. 0(c), we analyze the time T rev for the polarization 
of the central site to change sign. A strong dependence on 
the interaction range is found: T rev increases exponentially 
with L for a < 1, while for short-range interaction, it is in¬ 
dependent of system size. The exponential dependence of the 
reversal time on the system size is a consequence of the fact 
that long-range interactions induce a polynomially large en¬ 
ergy barrier, which must be overcome for the magnetization 
to reverse its sign. Such barrier induces ergodicity breaking in 
classical systems with long-range interactions (see Ref.[18] 
of the main text). On the quantum side, such barrier can be 



L 


FIG. 5: (Color online) Evolution of the polarization (a^(f)) for all 
sites n (top) and reversal time T rev of the central spin as a func¬ 
tion of the system size, L for different interaction ranges (bottom). 
Fized parameters: B = 1/2; W = 0, J z = 0. Initial state: 
Wx +1 )/ 2 (°)> = -1 and «*(l+ i)/ 2 (0)> = +1- Panels (a) and 
(b): a = 0; L = 9 for (a) and L = 13 for (b). Dark (black) curve: 
all sites n / (i + l)/2 and light (red) curve: n = (L + l)/2. The 
sign-reversal of the polarization occurs at the same time for all spins. 


overcome through the macroscopic quantum tunneling of the 
magnetization (see Ref.[17] of the main text). Since the en¬ 
ergy barrier is polynomial in the system size, an exponential 
tunneling rate follows, which explains the results presented in 
Fig. 0(c). Notice that T rev is distinct from the time T)/ 2 for 
the fidelity to reach the value 1/2. The latter grows polynomi¬ 
ally with system size. 


IV. LEAKAGE PROBABILITY 

Pieak refers to the probability for an initial state V), corre¬ 
sponding to an eigenstate of V, [see Eq. (1) in the main text] 
to leak outside its original energy band. The leakage proba¬ 
bility is defined as Pi ea k = 1 — lim.i_j.oo Pb(t), where Pb is 
the probability for the initial state to remain in the band b [see 
Eq. (3) in the main text]. 

If we start from a certain quantum mechanical state coupled 
with an amplitude e to another quantum mechanical state, the 
two being separated by an energy A, the probability to find 
the system on the second state will never be one; instead, for 
e/A -C 1, it will be at most of the order (e/A) 2 . Having this 
in mind we can estimate the asymptotic value of the leakage 
probability. The results for our estimates are valid in the dilute 
limit when the number of excitations is small compared to the 
system size, b/L <C 1. 


A. In the presence of an external field and no NN coupling 

For the Hamiltonian given by the Eq. (2) of the main text, 
with an external field and no NN coupling (J z = 0), each 
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state \V) in a band 6 is connected to approximatively L states 
in the nearest-neighboring bands. The coupling amplitude is 
h n , with (hi) = W 2 / 12. Thus, we can define the strength of 
the coupling as, 

W 

e ^7T T 

For small 6 ’s, the energy distance between the neighboring 
bands is proportional to JL (see main text), therefore 

r e 2 W 2 
Pleak OC ^ ~J2f j ' 

As an example, let us consider an initial state | V) from band 
6 = 1 , where all spins point up in the x-direction, except for 
one spin, which points down. The term of the Hamiltonian 
containing the external field can be written as, 


iff = J2 h n< = E + (6) 

n —1 n —1 

where the operators cr^' x flip the spins pointing in the x- 
direction. Due to this term, our initial state is connected with 
one state in band 6 = 0 and with L — 1 states in band 6 = 2, 
so that we have: 

P*=i~_fl + ( L _ 

Meafe — a 2 + 1+ 0,2 > 

^ 0,1 ^ 1,2 


where Ao,i = 2 J(L — 1) and A 12 = 2 J(L — 3) are the en¬ 
ergy differences between band 6=1 and the two neighboring 
bands. 

In general, starting with an initial state | V) in a band 6 , there 
are 6 connections with the 6—1 band and L — 6 connections 
with the 6+1 band, so that 


Pi 


leak 


~ 6 


A 2 

^ 6 - 1.6 


+ (L - 6) 


A 2 

M+i 


(7) 


where Ab-i^ = Eb_i—Eb = 2J(L— 26+1). This expression 
confirms the general scaling of Pi ea k given by Eq. 


B. In the presence of NN coupling and no external field 

We can estimate Pi ea k also for J z 7 ^ 0. The term of H 0 
containing the NN interaction can be written as, 

l- 1 l -1 T 

Hi- = E = E -i(°n' X +<' X ){°nZl+<£. i)- 

n=l n—1 

( 8 ) 

In general, J z connects a state in band 6 with m states inside 
that same band, n + states in band 6 + 2 , and n~ states in band 
6 — 2, such that m + n + + n~ ~ L. For a typical state with 
6 separated excitations, each one placed at least one site apart 
from the other, there are 26 states connected in the same band 
and + — 26 — 1 states connected in the outer bands, so that for 


b/L <C 1 the ratio n+ jrn oc L. Since the amplitude of the 
coupling is J z , we can write 

u n+ J? n_ J? J 2 
^-m + lA^/m+lAta/j^ 

Let us compute Pi ea k explicitly for the band 6=1. An 
initial state | V) from band 6 = 1 is coupled to two other states 
inside that band (apart from the situation where the initial state 
has an excitation on one of the border sites, 1 or L, in which 
case it couples only with one other state) and L — 3 states in 
band 6 = 3 (for an excitation on the border we have L — 2 
connections), so that we have: 

p b=i „ (L-3) J? (L - 3) J 2 Z 
leak ~ 2 A 2 a 8 J 2 (2L — 8) 2 ’ 


V. ENERGY SPREAD 


The time that it takes for the survival probability to reach 
the value 1/2 is denoted by Ti/ 2 - As discussed in the main 
text, the survival probability decay shows a Gaussian behav¬ 
ior, which justifies writing T \/2 ~ 1 /u>, where u> is the energy 
uncertainty of the initial state. Since the eigenstates of V in 
each band are degenerate, the perturbation Ho mixes them all. 
In this case, the energy uncertainty oj of the initial state can be 
approximated by the energy spread SE in each band induced 
by the perturbation. The analytical expression for the energy 
spread 5E evaluated at the first nonzero order of perturbation 
theory is studied below. 

We compute the eigenvalues of the total Hamiltonian given 
by Eq. (2) of the main text, using second order perturba¬ 
tion theory for degenerate levels. Note that our unperturbed 
Hamiltonian is the long-range part V, while we consider 11 0 
as the perturbation. In the following, we set £>, J z = 0, so that 
the perturbation Hq is determined only by the random field, 
as in Eq. ([ 6 ]) above. 

Let us consider the energy band 6=1. We denote the initial 
state \V) in this band as |1, k ), where k = 1,..., L indicates 
the position of the excitation. In labeling the states, we neglect 
the double degeneracy due to the flipping of all spins, since the 
states with M x = (L — 2)/2 are only connected to those with 
M x = — (L — 2)/2 in a very high order of perturbation theory. 

The perturbation [Eq. @ above] connects the initial 
state in band 6 = 1 to the state |0, 0), belonging to band 6 = 0, 
and to L — 1 states in band 6 = 2. The latter states are denoted 
by |2, k,j)\ they have one excitation on site k = 1,..., L — 1 
and the other on site j = k + 1,..., L, so that the total number 
of states is L(L — l)/2. 

Since the degeneracy inside the band is not removed at first 
order of perturbation theory, we use second order perturbation 
theory for degenerate levels, namely the eigenvalue problem, 

(V + eH™)\4> 1 } = (Ei + eE{ + e 2 ^)!^) (10) 
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where 

L 

IVh) =^c o , fc |l,fc) + (ec^, o + e 2 c^ 0 )|0,0) (11) 

k=1 

L—l L 

+ E E ( ec +,fcj + e2c +^fc,j)|2, k,j). 

k =1 j'=fc+l 

The action of the “unperturbed” Hamiltonian on the unper¬ 
turbed states is trivial, 

V\b,k) = E b \b,k), with 6 = 0,1,2,... 

while the perturbation acts as, 

K\k k) = h k |0,0) + E hj\2, k,j). (12) 

&k 

Collecting the e terms in Eq. (flOt gives, 

L 

-E[ c o,fc11 7 k) + (Eq - £i)d, o |0,0) + 

k=l 

L — l L 

(E 2 -E 1 )J2 E c Uj\ 2,k,j) + 

k=l j—k-\-l 

E c o,fc f^| 0 , 0 )+E^| 2 ,^i) ) =0 (13) 

fc=i y ) 

Bracketing Eq. ( fT3] > respectively with (0,0|, (1, s | and 
(2, a, /3\ with /3 > a, we obtain. 


V ’° E 1 - E 0 


E[ = 0 , 


E c o,khk , 


(14) 




C +,a,/3 — tt, tt. ( c 0 ,ahp + CQ t ph a ) . (15) 

-C/1 — £/2 

Collecting the e 2 terms in Eq. ( I I ()[ i and taking into account that 
E[ = 0, we have, 

L 

(Eq - £i)d o |0,0) - El 1 E c o,fc|l, k) + 

k=l 

L — l L 

(e 2 -e 1 )J 2 E 4 EI 2 -M + 

fc=l ji —/c+1 

L—l L 

c 1 -, 0 ^ 10 , 0 )+ E E c +,k,jH(^ 1 2 , fc, j) = 0 . 

/c=l J = /c+l 

(16) 

Bracketing Eq. ( IT6t with (0,0|, (1, s| and (2, a, /3| with /3 > 
a, we get d 0 = 0, a p = 0, and 

E{ J c o, a = ci i0 (l, |0,0) + E S I#<T |2, fc, j), 

j>k 


The equation above, due to the symmetry of the coefficients 
c+ k j = c r + J k [see Eq. ( flsT i ] can be rewritten as. 


E^co^s = cLfihs + E c +,s,j h r ( 18 ) 


Inserting Eq. ( im and Eq. ( fT5l) in Eq. (fl8l ). one finds that the 
second order corrections E( 2 are the solutions of the L equa¬ 
tions, 


CQ,s 


E 




(Ei - E 0 ) “ (Ei - E 2 ) 


- E 


ii 


+ 


El cojhghj 


(Ei — Eq) (Ei — E 2 )_ 


= 0 , 


(19) 


with s = 1, L. In other words, E( 1 are the eigenvalues of 
the symmetric matrix C, whose diagonal elements are. 


Css = 


h 2 


E- 


hi 


(El - Eq) ^ (Ei - E 2 ) ’ 


( 20 ) 


and off-diagonal elements for k ^ s are, 

1 1 


Cks — hkh s 


(Ei — Eq) (Ei — E 2 ) _ 


( 21 ) 


Let us now estimate the eigenvalues of this matrix in the 
limit of large L, at fixed IT'. We know that, 

Ei - Eq = 2 J(1 - L) and Ei - E 2 = 2 J(L - 3), 

therefore, in the limit of large L one has that the off-diagonal 
elements are negligible with respect to the diagonal ones, 
since Cks ~ o(l/L) for k ^ s, while C ss ~ o(l). 

Thus, we can estimate the eigenvalues from the diagonal 
elements only. In particular, since we are interested in the 
energy spreading, we evaluate 


6E' = (Ci s ) w - (C ss )l 

W 4 ( 1 L—l 


> 2 
! w 


180J 2 \ (L- l) 2 (L-3) : 


( 22 ) 


where in computing the average over disorder we took into 
account that 


( h s h k ) w =0 for s ^ fc, 


and 


, 2 , - — lh 4 \ - — 

\' l S/W 12 ’ \"'S/W gQ • 

In the limit of large system size, we therefore get the 
asymptotic behavior, 

W 2 


jVl' 


(17) 


SE 


(23) 
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The generalization to an arbitrary band b is far from triv- which gives the same estimate as in Eq. (|23T >. 
ial. We will provide the details of this derivation in a longer 
version of this paper. We present here only the final result. 

Similar to the case b = 1, we can estimate the energy spread¬ 
ing for the general b band as, 






